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ABSTRACT 

We present Chandra X-ray monitoring of the M87 jet in 2002, which shows that the intensity of HST-1, 
an optical knot 0.8" from the core, increased by a factor of two in 116 days and a factor of four in 2 yrs. 
There was also a significant flux decrease over two months, with suggestive evidence for a softening 
of the spectrum. From this variability behavior, we argue that the bulk of the X-ray emission of 
HST-1 comes from synchrotron emission. None of the other conceivable emission processes can match 
the range of observed characteristics. By estimating synchrotron model parameters for various bulk 
relativistic velocities, we demonstrate that a model with a Doppler factor 5 in the range 2 to 5 fits 
our preliminary estimates of light travel time and synchrotron loss timescales. 

Subject headings: galaxies: active — galaxies:individual(M87) — galaxies:jets — magnetic fields — 
radiation mechanisms:nonthermal 



1. INTRODUCTION 

With the advent of Chandra X-ray observations of ra- 
dio jets, evidence supporting the synchrotron process for 
the origin of the X-ray emission from the jets of low lu- 
minosity radio galaxies has grown. One line of evidence 
is the similarity of the X-ray morphology to the radio 
and optical structures which are most likely synchrotron 
emission because of the detection of linear polarization. 
Another argument is that the values of the spectral in- 
dex of power law fits to the X-ray spectra are generally 
steeper than the radio spectra whereas inverse Compton 
(IC) models predict ax < dr'^ (e.g. Hardcastle et al. 
2002; Harris & Krawczynski, 2002). 

In this paper we exploit another fundamental differ- 
ence between IC and synchrotron emission: the variabil- 
ity time scale. Although the injection of new relativistic 
electrons can produce the same rise time for both pro- 
cesses, the characteristic loss timescales are vastly dif- 
ferent. All IC models require relatively low energy rel- 
ativistic electrons. For synchrotron self Compton (SSC) 
models, Lorentz energy factors, 7, are normally compa- 
rable to the values of a few thousand which are respon- 
sible for the observed radio emission. For IC scattering 
between the relativistic electrons and photons of the cos- 
mic microwave background (IC/CMB), 7 <1000 (Harris 
& Krawczynski 2002). Synchrotron models however re- 
quire 7 K, 10*" with corresponding half-lives of order a 
year for the sorts of magnetic field strengths commonly 
estimated for features in radio jets. 

Our data show strong intensity variability from the 
unresolved core and HST-1. We use the flux doubling 
time to derive an upper limit to the source size and the 
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decay time to constrain the timescale for E^ losses, and 
in turn, 5. We take the distance to M87 to be 16 Mpc so 
that one arcsec = 77pc. 

2. THE X-RAY MONITORING OBSERVATIONS 

In 2002 we obtained five observations of 5ks each on 
Janl6, Febl2, MarSO, JunOS, and Jul24. To minimize 
pileup, we used a readout time of 0.4s for the stan- 
dard l/8th subarray on the ACIS-S3 chip (PUG 2001). 
Data reduction followed the 'threads'^ (CIA02.2). To 
recover the inherent resolution of Chandra, we removed 
pixel randomization and rebinned the data to 0.1 native 
ACIS pixels. For 3 energy bands (soft, S=0.2-0.75keV, 
medium, M=0.75-2keV, and hard, H=2-6keV) flux maps 
were obtained by standard procedures; a correction was 
applied for degradation of the quantum efficiency by 
measurement of fixed regions of hot gas emission near 
the jet. A wide band image was made by adding the 3 
energy bands and fig. ^ shows the result of averaging our 
fluxmaps. 

Our monitoring shows striking variability in both HST- 
1 and the unresolved core (which could well come from 
jet segments closer than 0.4" to the nucleus). The other 
knots show weak or no variability and will be analyzed in 
a future paper. Contour diagrams of the core and HST- 
1 for each observation can be found in Harris (2003a); 
here we show lightcurves for the core and HST-1 (fig. Ell. 
The photometry was obtained by taking small circles 
with radii 0.44" and subtracting the background from 
rectangles above and below the jet (with a total area of 
13.6 arcsec^). 

The largest increase occurs for HST-1 between Mar and 
Jul with over a factor of two brightening in 116d^. This 
implies that the emitting volume of the variable compo- 
nent has a characteristic size < the light travel time^*' 
which is O.lpc for a stationary source. The brightening 
in the core between our last two observations showed a 
similar rate of increase with a 20% gain in 46 days. 

* http:/ /asc. harvard.edu/ciao/threads/ 

^ A brightening was also seen in our HST monitoring between 
2001 July and 2002 July. Those data will be analyzed in a separate 
paper (Biretta et al. 2003, in preparation) 

Although see Protheroe (2002) for caveats. 
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Fig. 1. — Chandra X-ray image of the M87 jet constructed by averaging the data of Wilson &: Yang (2002) taken in 2000 July, together 
with our 5 observations. The contours increase by factors of two in brightness, with the lowest contour level being Ix 10~^^ erg cm~-^ per 
pixel in the 0.2 to 6 keV band. Jet features are labelled according to the standard convention and a 10" scalebar is shown. ACIS pixel 
randomization has been removed; the data were regridded to a pixel size of 0.0492"; and a Gaussian smoothing function of FWHM=0.25" 
was applied. 
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Fig. 2. — Lightcurves for the core(left panel) and HST-1 (right panel). The curves are for different energy bands: total (solid line), hard 
(broken), medium (dotted), and soft (dashed). The data at the left edge come from the archival 2000 July observation (Wilson and Yang 
2002) and those to the right from our 2002 monitoring. 



The relative increases in the three energy bands 
demonstrate a hardening of the spectrum for HST-1 be- 
tween March and July 2002 (compare the rates of in- 
crease in each band in fig.|21). The intensity increase was 
a factor of 2.7±0.3 for the hard, 1.7±0.1 for the medium, 
and 1.6±0.1 for the soft band. Although we defer spec- 
tral analysis to a later paper, this hardening produced a 
change in the X-ray spectral index from ax ~ 1-35 ± 0.1 
(prior to 2002 May) to = 1.05 ± 0.07 (for 2002 July). 

The 'cleanest' drop occurs for HST-1 between the 2002 
January and February observations. The intensity ratios 
(Feb over Jan) for these 27 days are: 0.79±0.09 (H), 
0.83±0.04 (M), and 0.87±0.05 (S). While these values 
are the same within the errors, there is a progression con- 
sistent with expectations for losses (but contrary to 
that of increased absorption): the hard band is dropping 
faster than the soft band. A similar effect is seen in the 
core decay between Mar and Jun: the hard band dropped 



32±7% while the medium band dropped 14±5%. 

All of these timescales are subject to considerable er- 
ror. Although the statistical errors for the flux values 
are less than 10%, there will always be some flux present 
which is not part of the variable component and our sam- 
pling interval is not small enough to determine shorter 
variations. 

Variability of a synchrotron source can result from a 
number of effects such as changes in any of the absorption 
along the line of sight, the magnetic field strength, B, 
the number of particles contributing to the radiation in 
a given band, or the beaming factor (Protheroe 2002). 

Adiabatic compression or expansion is particularly ef- 
fective in changing the radiated power since the electrons 
gain or lose energy and at the same time, B changes to 
augment the effect. Given a power law electron energy 
distribution N(E)(xE~p, there should be no change in p. 
On the other hand, if the effects of a high energy cutoff 
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extend down to the energy range of interest, a compres- 
sion would cause an increase of intensity accompanied 
by a spectral hardening since the emission from a fixed 
observing band would come from lower energy electrons 
and the spectral curvature would be shifted to higher en- 
ergies. However, for the remainder of this paper we focus 
on an alternative synchrotron model based on changes to 
the number of radiating particles. 

For this model, increased intensity arises from the in- 
jection of new radiating particles with a spectrum flatter 
than that which existed previously. The spectral soften- 
ing which accompanies flux decline comes from normal 
losses. The chief difference in these models is that 
the compression/expansion model requires that N(E) de- 
parts from a power law. Our data allow spectral curva- 
ture but the statistics are not good enough to demand 
curvature. 

3. THE SPECTRAL ENERGY DISTRIBUTION OF HST-1 

In order to derive synchrotron model parameters, we 
need to know the overall spectral distribution of the emit- 
ted energy. Both the radio and optical morphologies 
of the inner jet are shown in fig. 3 of Perlman et al. 
(1999) with 0.2" resolution. At higher resolution (fig. 
1 of Biretta, Zhou, & Owen 1995) it can be seen that 
there are a series of several knots between 0.8" and 1.2" 
from the core. HST data demonstrate that the bright- 
est upstream knot is slowly moving but the downstream 
knots are moving with an apparent speed of 6c (Biretta, 
Sparks, & Macchetto 1999). Our feature 'HST-1' (fig.IH 
corresponds to the upstream, optically brighter part of 
this segment of the jet. We take the radio flux density 
corresponding to the X-ray knot to be 3.8 mJy, as mea- 
sured on the high resolution 15 GHz map. From a new 
VLA observation,^^ we derive an upper limit of 2 mJy at 
43 GHz. 

The optical flux densities come from the detailed spec- 
tral analysis at 0.15" resolution based on data obtained 
in 1998 February by Perlman et al. (2001). The spec- 
trum is shown in fig. 01 and for the purposes of sec.0I we 
approximate these data with a single power law between 
10^ and 10^^ Hz with q;=0.68 and amplitude determined 
by S(2keV)=100nJy, parameters which describe the ob- 
servables to within a factor of two. 

4. A 'modest beaming' synchrotron model for hst-1 

Previous analyses have demonstrated that the X-ray 
jet emission of M87 is not thermal emission (Wilson & 
Yang 2002; MarshaU et al. 2002; Biretta, Stern, & Harris 
1991); is not SSC emission (ibid.); and is not IC/CMB 
emission because of the large values of F and small angles 
to the line of sight required (Harris & Krawczynski 2002). 
The X-ray spectra of the knots are much steeper than in 
the radio or optical and strongly favor synchrotron radi- 
ation (Wilson & Yang 2002), confirming the conclusion 
of Biretta, Stern, & Harris (1991). Our variability data 
is fully consistent with X-ray synchrotron emission and it 
remains only to demonstrate that plausible physical pa- 
rameters can be found with a simple model which posits 
that the jet has a relativistic bulk velocity, as implied by 

We obtained a short observation in 2002 October as a result 
of an ad hoc proposal. The VLA was in C array with a beamsize 
of 0.5". 
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Fig. 3. — The spectrum of HST-1. Optical flux densities come 
from HST data (1998) of Perlman et al. (2001), but with smaller 
apertures to correspond to just the first, brighter, part of the com- 
plex. The X-ray flux densities come from the observed fluxes (i.e. 
no correction for absorption has been applied) assuming ax=l. 
Values are shown for the faintest (2000 July) and brightest (2002 
July) intensities. The actual flux densities may be somewhat larger 
since we used a small aperture for photometry to isolate the vari- 
able fraction of the total intensity. 



the one-sided emission at all frequencies. In this section, 
wc do that for HST-1. 

Using the time required for the flux to double pro- 
vides an upper limit to the source size. To use the 
observed decline, we take the observed variables back 
to the jet frame, calculate the properties of the syn- 
chrotron source assuming equipartition conditions (Pa- 
cholczyk 1970), and then evaluate the halflife of the elec- 
trons responsible for the X-ray emission, both in the jet 
frame and as observed at the Earth. We also compare 
energy densities in the magnetic field and in the syn- 
chrotron photons. We performed these calculations for 
8 values of S between 1 and 16. For converting various 
quantities to the jet frame and back, we used expressions 
from the appendix of Harris and Krawczynski (2002). 
The results for representative values of S are shown in 
Table HI 

We have made model specific estimates of how to relate 
an observed decline of 21% (in the hard band in 27 days) 
to the effects of E^ losses on the electron distribution 
(Harris 2003b). We find consistency for S in the range 
2 to 5. However, the main conclusion is that regardless 
of the precise model and the actual values found for the 
rate of decline, the observed timescales for decay arc fully 
consistent with the beaming parameters for a jet angle 
of order 10° to 17° (to the hne of sight) suggested by 
the HST proper motion studies of Biretta et al. (1999) 
in their Table 3. Larger beaming factors imply weaker 
fields and longer halflives, and vice-versa. A recent HST 
observation shows the varying region to be unresolved in 
the optical/UV, independently placing an upper limit of 
0.02" (1.5pc) on its size (cf. Table[Il). 

Values of 5 in the range 2 to 5 involve very modest 
synchrotron luminosities in the jet frame, equipartition 
magnetic field strengths of order 1 mG, and a total en- 
ergy stored in particles and fields which can be supplied 
by a low drain on the jet kinetic energy, « 10^" erg s~^, a 
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Table 1. Beamed Synchrotron Parameters for HST-1 



Parameter 


<5=1 


S=2 


<5=4 


<5=6 


<5=8 


(5=16 



radius (arcsec) 
K (erg s-l) 
log Kot (ergs) 

W (mG) 

u'(B) (erg cm~'^) 
u'(sync) (erg cm~'^ 
r; (yr) 

Ti (Earth) (yr) 



0.0012 


0.0024 


0.0048 


0.0072 


0.0096 


0.0192 


6.5E40 


4.1E39 


2.5E38 


5.3E37 


1.6E37 


1.0E36 


48.55 


48.40 


48.17 


48.04 


47.96 


47.75 


135 


181 


244 


274 


328 


443 


4.2E6 


5.7E6 


7.7E6 


8.7E6 


1.0E7 


1.4E7 


13 


3.7 


1.0 


0.51 


0.28 


0.077 


6.9E-6 


5.4E-7 


4.0E-8 


l.OE-8 


3.7E-9 


2.4E-10 


7.5E-6 


1.2E-7 


1.8E-9 


1.9E-10 


2.9E-11 


4.5E-13 


0.011 


0.126 


1.3 


4.3 


12.1 


81.9 


0.011 


0.063 


0.32 


0.72 


1.5 


5.1 



Note. — The radius is the angular size corresponding to the light travel time in the jet frame. Jet frame parameters are primed. L^, 
is the synchrotron luminosity integrated over the band which is observed at the Earth between E9 and E18Hz. Etot is the total energy 
in particles and fields required to produce L^. 71 2 are the lower and upper bounds on the electron energy distribution. The u's are the 
energy densities in the magnetic field and the synchrotron photons, ti is the halfiife of electrons producing the highest energy photons. 



value <1% of the jet power estimated by Young, Wilson, 
& Mundell (2002). With these parameters, the effective 
energy density of the CMB as seen by the jet is less than 
the synchrotron energy density for most choices of al- 
lowed r, and the combined photon energy density <0.2 
of that in the magnetic field for all (5 >2. Thus the bulk of 
the energy loss is via the synchrotron, not the IC channel. 
The model predicts that the optical emission will decay 
with a timescale of order ten times longer than that we 



observe for the X-rays since T| oc ^. For a straight jet, 
an angle to the line of sight oi 6 ~ 10° implies a total jet 
length of just under 10 kpc. 

These estimates are based on decay timescales and our 
results would change if the actual decay time is shorter 
than our current estimate, in which case we would choose 
a smaller value of 5; if the magnetic field is substantially 
less than the equipartition value and thus ri would be 
larger than calculated because of the weaker field and we 
would again choose a smaller 5\ or if the source size is 
substantially less than the light travel time which would 
mean that the field is larger than calculated, halfiives 
would be shorter, and we would select a larger 5. 

5. CONCLUSIONS 

We find that a synchrotron model with a beaming fac- 
tor of (5=2-5 gives a good fit to our observations for 
HST-1 and is compatible with radio and optical inter- 
pretations. It remains to be seen if similar models can 



be constructed for powerful jets such as that in 3C 273. 
Many of the gross features such as the progressions of 
X-ray emission being stronger close to the core, whereas 
the radio intensity increases moving out along the jet; 
and the slight shifts in brightness (radio brighter down- 
stream of the X-ray peak for some knots) are the same 
in M87 and 3C273. It is also worth noting that if M87 
were further away or observed with poor resolution, the 
variability which we now know arises from a knot in the 
jet might well be ascribed to conditions very close to the 
black hole (as in blazar variability studies). We shall con- 
tinue to monitor the jet of M87 with Chandra between 
2002 November and 2003 July with 8 observations and 
the jet will also be monitored at similar intervals with 
the HST. These data should provide refined estimates of 
variability timescales, spectral evolution for the core and 
HST-1, and variability in the other knots. 
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